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We derive the probability for rapidity gaps in a parton cascade and investigate the 
dual connection with hadronic final states. A good description of observations in 
e+e~-annihilations is obtained by perturbative QCD calculations in MLLA using 
previously determined parameters (QCD scale A and fcy-cut-off Qo) and applying 
the parton hadron duality picture. Further predictions are derived; especially, for 
gaps between jets at variable resolution we predict a strong variation of gap proba- 
bilities for small parameters ycut — > in the transition from jets to hadrons. Large 
gaps between partons correspond to large spatial separations of colour charges: a 
colour blanching mechanism by soft processes is suggested. 



1 Introduction 

The occurrence of large rapidity gaps in e"*'e~-annihilations, or more generally, 
inside quark and gluon jets, provides an interesting testing ground for models 
of colour confinement and hadronization. 

In case of e^e~-annihilation the primary process is the production of a 
qq pair. According to a perturbative mechanism aa large rapidity gaps in 
the hadronic final state occur if in a subsequent process two low mass parton 
pairs arc formed in colour singlet states {qq or gg) recoiling against each 
other. The production rates predicted for these "hard" colour singlets which 
involve highly virtual intermediate gluons, however, are much smaller than 
those observed by SLDJS by about two orders of magnitude. ^ 

A quantitative description of the data is provided by the JETSET MC □ 
which combines an initial parton cascade, cut off at a scale ~ 1 GeV, with 
a string hadronization model. The SLD data are shown in Fig. la. ^ The 
JETSET result fits the data in the full range of rapidity intervals Ay B; for 
large Ay the contribution of r+r" events becomes important and the dashed 
histogram represents the purely hadronic contribution. 

We will discuss here another approach, originally proposed for the treat- 
ment of single inclusive spectra ("Local Parton Hadron Duality"- LPHEB) 
in which the parton cascade is perturbatively evolved further with cutoff 
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Figure 1. Probability for raiw|dity gaps of length l\y\ (a) between charged particles in e+e~ 
annihilation: i^ata from SLDcl, also shown is the JETSET MC for qq initial state without 
t'^t~ events.tl Furthermore we show our results from ARIADNE MC at the parton level 
which, in the duality picture, correspond to gaps between all hadrons (curves I, II); an 
estimate for gaps between charged particles is obtained by multiplying these results with 
the ratio of gap fractions f^h ^ jrall fj-Qjjj ^j^g -j^q after string hadronization (curves 
III,IV) (Parameters used in MC: (I) A = fl.20 GeV, A = 0.015, Nj^ = 3; (II) A = 0.32 
GeV, A = 0.015, = 5). The data pointsH are moved to the right edge of every interval 
as appropriate for a cumulative quantity; (b) vertical bars represent the same ARIADNE 
results (curves I, II in (a)) which are also close to the data; the curves (I) here correspond to 
our analytical calculation of the Sudakov form factor, the curves (II) include also the cross 
over correction (full curves: A = 500 MeV, Durham fcy; dashed: A = 350 MeV, standard 
fcy; A = 0.015 always) 



kr > Qo Sit Qo > A with QCD scale A ~ few 100 MeV. The resuhs are com- 
pared directly to the data without any explicit hadronization phase and the 
partonic final state is assumed to represent the hadronic one in a dual sense. 
This simple approach has been applied to a variety ^nf problems, mainly on 
inclusive quantities, with rather surprising successes .ltlI 

Our present study carries this idea further, in that we consider an ob- 
servable which tends to become exclusive in the limit of large gaps. Indeed, if 
colour charges are separated by a large gap, then one might expect, accord- 
ing to conventional wisdom, a neutralization by non-perturbative processes 
leading to large deviations from the perturbative calculations. At the time of 
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a previous study along these lineaJ the ratio Qo/A was not weU determined 
and only an upper limit of the gap fraction has been_given. Meanwhile, from 
an improved analysis of jet and hadron multiplicitiesu a determination of this 
ratio has been obtained. Now, the gap rate can be predicted in absolute terms 
within the perturbatively based duality picture. 

2 Rapidity Gaps in the Perturbative Parton Cascade 

The probability for no radiation into a certain angular interval is given in 
field theory by the exponential Sudakov form factor originally derived in 
QED. In our application we consider the rapidity gap without gluons above 
the transverse momentum cutoff Qq. Let us consider specifically the angular 
interval between Oi and 82 (Oi > 62), the rapidity is then obtained from 
y — — Intg^. Let us further denote the probability for emission of a gluon 
at an angle 0' with the energy lo' off a parent parton p (either a gluon(g) or 
a quark(q)) as pp{Lj',Q') = dup / duj' dQ' . Then, the Sudakov form factor for 
the angular ordered cascade is given by 

ApiP, e, Qo) = exp(-u;p(P, 9, Qo)) (1) 

WpiP,&,Qo) = [ duj' [ de'ppiu;',e'), (2) 

J Jk±>Qa 

and it represents the probability for no gluon being emitted within the cone 
of half angle at transverse momentum above Qo from the parent parton p 
of energy P = Q/2. In particular, 

Ap{e2)/Ap{ei) - exp (-uv(02) + M^i)) (3) 

represents the probabilty that there is no emission of a gluon with emission 
angle between 8i and 82. These rates have been calculated in different 
approximations 

The double logarithmic approximation (DLA) 

The simplest approximation takes into account only the leading contributions 
from the angle and energy singularities of the gluon emission. The gap proba- 
bility fp is easily calculated analytically! For the symmetrical gap in the cms 
we find a good approximation for not too large gaps Ay/2 ^ Y 

fp{Ay) ~ exp(-ApAy) (4) 
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with Cg = 3, Cq — ^. The gap rate decreases exponentially with Ay. The 

slope depends sensitively on A for small A and A — > oo for A — > 0. 

The modified leading logarithmic approximation (MLLA) 

In this improved approximation also the next to leading loaarithmic terms are 

included. Some analytic results have been obtained before.! Here we calculate 

the probability Wq as in the multiplicity analysis □ 

Wq = dn' dz^^^q,{z) (6) 

by numerical integration, where k — Q sin(0/2) is the jet virtuality at opening 
angle 0, the sphtting function is ^qg{z) = 2Cf(1 + (1 — 2)^)/z and fcr denotes 
the transverse momentum; this is taken as fcr = z(\ — z)k ("standard") or 
fcy — niin(z, 1 — z)k ("Durham"). The different kr lead to somewhat different 
A without changing A = 0.015. 

In this calculation the exponent Wq is of 0{as). To next order, processes 
play a role where a secondary gluon is emitted into the gap although the firstly 
emitted gluon is outside the gap, this we call "cross-over effect" . We note that 
there is no such effect for 8' < 82/2 because of angular ordering, therefore, 
the maximal effect is a shift by Ay = In 2 '--^ 0.7. Otherwise, we take the effect 
into account by a correction factor in DLA accuracy: we multiply (H) under 
the integral with the probability P for no secondary emission into the gap, 
approximately with 

P=exp(-i(u;,(fc',8')-^«,(fc',82-8'))) for ^ < 8' < ^ (7) 

P=exp(-^{wg{k',Qi^e')-Wg{k',e2-Q'))) for & > ^ (^) 
The Parton Monte Carlo 

As axontrol of our analytical calculations we compare also with the ARIADNE 
MCLJ at parton level which is based on similar principles, i.e. cutoff kr > 
Qo and possibility to choose a small A parameter. We take parameters (I) 
determined from a fit to hadron multiplicitiesEj; new parameters (II) are 
determined to improve the fit to jet multiplicities at small ijcut- In the MC 
we used only it-quarks with niu — and kept number of flavours Nf fixed. 

3 Comparison with Data 

The MC results for the two sets of parameters (I) and (II) are also shown in 
Fig. la. They refer to the gaps between all hadrons in the duality picture. 
To obtain an estimate for the gaps between charged hadrons only we multiply 
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Figure 2. Probability for rapidity gaps of size h.y between jets for different resolution 
parameters ycut = QT^^ /s in the Durham algorithm, so Q%'^^ = corresponds to full 
resolution, i.e. the hadronic final state. The full line represents ARIADNE MC (curve II 
in Fig. 1), the others the analytical Sudakov calculations without crossover effect (dashed: 
A = 350 MeV, standard kr; dash-point: A = 500 MeV, Durham fey, A = 0.015 always, 
Q=91 GeV) 



these curves with the respective ratio derived from the parton and full hadron 
MC which we parametrized as f^^ j jaii ^ i _^ 2 Ay - 0.1(Ay)2. After this 
correction, one observes a very good agreement of this 2-parameter model 
with the data (at large Ay one should compare to the dashed line). 

In Fig. lb these MC results are represented again as vertical bars for 
reference, also to the data. The curves (Il-represent our analytical calculations 
based on (^) for two sets of parameteraa, curves (II) include also the cross 
over corrections (|7|),(^). Good agreement of the latter results with the MC is 
obtained up to Ay '--^ 3 — 4, it falls below the MC at higher Ay but could still 
be close to the experimental data. 

4 Further Predictions 

We note three consequences of our approach which follow directly from the 
simple DLA formula (^). 
Quark vs. gluon jet 

The slope is proportional to Cp as the Sudakov form factor is derived from 
the 0{as) gluon emission probability. Then, in DLA, the slope in a gluon jet 
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is larger by Ca/Cf = 9/4 as compared to a quark jet. 
Energy dependence 

The slope behaves like A ~ lnln(P/(5o) so the gap distribution gets steeper 

with increasing jet energy P. 

Dependence on cutoff Qq and jet resolution 

In the duality picture, the cutoff Qo appears as a hadronization scale which 
limits the resolution of separate partons. In the evolution equation it can be 
interpreted also as jet resolution parameter in the Durham algorithm whith 
Ucut = (Qcut/Q)'^ with Qo replaced by Qcut- Therefore, we expect a strong 
dependence of the slope on the jet resolution parameter through A ^ — In A. 
In the theoretical calculation all hadrons are resolved for Qcut Qo in the 
duality picture, experimentally for Qcut ~^ 0. This mismatch can be resolved 
by relating^ {Q*^tY = {QlulY + Ql- In Fig. 2 we show predictions for the 
rapidity gap probability refering now to jets at resolution Q^Mt" = ^JVcutQ in 
the Durham algorithm; we neglect the crossover effects which should be small 
at larger Qcut- One can see the dramatic rise of f{Ay) at Ay = 4 by about 
three orders of magnitudes if we replace hadrons by jets at resolution 1 GeV. 
It will be interesting to verify this new effect. 

5 Conclusions, a Puzzle and a Physical Picture 

We have derived perturbative predictions for rapidity gap distributions using 
the two parameters Qq and A from earlier fits to the mean global particle 
multiplicity. The agreement with the SLD data is quite remarkable and in 
support of the simple duality picture also in case of this new, partially exclu- 
sive, observable. It would be desirable to determine gap fractions from final 
states with inclusion of neutral particles which would allow a more direct 
comparison with our calculations. A crucial test of our picture is the strong 
dependence of the gap distribution on the jet resolution. 

Whereas the phenomenological description of the model is successful, the 
interpretation imposes a serious puzzle: for a large gap, take Ay — 3, the 
first gluon in each hemisphere is emitted only after a mean lifetime of about 
10-20 f, so the perturbative evolution is not disturbed, even if the initial qq 
pair gets separated far beyond the typical confinement distance]^ 

As the mechanism for global colour blanching is not quantitatively known 
we may consider phenomenological scenarios. The colour blanching could be 
mediated by "gluers"l3 at kx '-^ Qo, they are expected to cause the produc- 



" We estimate the lifetime of the virtual quark radiating a gluon with momenta k, kx as| 
r ~ — — — ~ i and take the average of r over ui' and 0' as in (bl) within the DLA. 

771 g tTlq Krj-i Lj 
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tion of hadrons with flat rapidity plateau, already in absence of perturbative 
gluons. In our interpretation there should not be any associated real hadron 
production in the blanching process as the gap would be refilled. The con- 
finement effects must be weak enough so that the successful perturbative 
calculation of the gap rate is not invalidated. We consider the possibility that 
in the field of the separating qq pair (say, inside a tube with virtual partons 
of kx < Qa) a very soft new qq pair is produced to ensure confinement at 
usual distances, and further on at other vertices (Fig. |^). In a more spe- 
cializedptnodel the picture in Fig. ^ could be realized by effective hadronic 
verticesll^. At any rate, it will be interesting to test further the predictions of 
the perturbative analysis at low energy scales. 




Figure 3. Hadron final state emerging from qq with large rapidity gap: the first gluon 
emission occurs far outside the confinement region of size ~1 f. A possible mechanism is 
colour blanching by soft qq pairs with fcy ^ Qo at all vertices. 
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